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Abstract 

We investigate the structure, chemical bonding, electronic properties, and magnetic behavior of 
a three-dimensional graphitic network in aba and aaa stacking with intercalated transition metal 
atoms (Mn, Fe, Co, Ni, and Cu). Using density functional theory, we find induced spin-polarization 
of the C atoms both when the graphene sheets are aba stacked (forming graphite) and aaa stacked 
(resembling bi-layer graphene). The magnetic moment induced by Mn, Fe, and Co turns out to 
vary from 1.38 {ib to 4.10 whereas intercalation of Ni and Cu does not lead to a magnetic 
state. The selective induction of spin-polarization can be utilized in spintronic and nanoelectronic 
applications. 
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I. INTRODUCTION 



The rich variety of the physical and chemical properties of nanoscale materials is used in 
a wide range of advanced technological applications. Graphene, a single layer of graphite 
with C atoms tightly packed in a honeycomb lattice is one of the most promising low- 
dimensional nanomaterials, attracting immense interest of both experimentalists and theo- 
reticians because of its two-dimensional (2D) structure and unique properties. In both 
pristine graphene and graphite no magnetic ordering is expected, by experiment and theory. 
However, there are various experimental methods to induce magnetism: hybridization be- 
tween the C p z and Ni d states [4(1, ion implantation [5], proton irradiation of highly oriented 
pyrolytic graphite jf], Q], and point defects jsj]. Ferromagnetism can be induced in multilay- 
ered graphene by randomly removing single C atoms [9] . Theoretically, spin-polarization has 



been confirmed for various C defects 



1CH16I] , adsorption of molecular oxygen and hydrogen 



on graphene |l7|] , and substitutional Mn doping of graphene [18| . 

Inducing spin-polarization in graphene by doping transition metal (TM) atoms [19M21] 
is important because it can lead to scattering [22| and modify the electronic states locally, 
which is required for graphene-based electronics and Kondo physics. Often, very simplistic 
assumptions about the effects of the TM atoms are not valid, but the properties depend 
strongly on the host [23]. Furthermore, precise knowledge of the TM-C interaction is impor- 
tant for understanding carbon nanotube growth [24] , fuel cells [25] , and the role of implanted 
magnetic atoms, such as Fe, for the magnetic order j^fj. For these reasons it is surprising 
that TM intercalation in a three-dimensional (3D) graphitic network has not been studied 
with respect to the modifications of the interlayer interaction. In this paper, we will present 
first-principles results on the spin-polarization induced by Mn, Fe, Co, Ni, and Cu atoms 
intercalated between adjacent C layers. We study two types of stacking: aaa, which results 
in multilayer graphene, and aba, which results in bulk graphite. The different d valences 
of the TM atoms result in a range of induced magnetic moments (magnitude and spatial 
distribution) as needed for spintronic applications. 
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II. METHODS 

Calculations based on density functional theory [27] are carried out using the generalized 
gradient approximation (GGA) and the plane wave Quantum-ESPRESSO code [28|], which 
already has been used successfully for describing graphene A strong corrugation 

of the crystal potential is present in layered materials in the direction perpendicular to the 
atomic planes. Consequently, the local density approximation is not valid. Even though the 
van-der-Waals forces are neglected in both the local density approximation and the GGA, 
the GGA leads to a reasonable matching of the interlayer distances with experiments and 
is adequate for systems with an inhomogeneous charge density We apply Perdew- 



Burke-Ernzerhof pseudopotentials [38] and include an onsite Coulomb interaction correction. 
In our calculations scalar relativist ic effects are included, while spin-orbit coupling is not 
taken into account as it is negligible in the present systems. The value of the onsite interac- 
tion U is varied between and 8 eV. Our results show that it is essential to use a non-zero U 
for the d electrons of the TM atoms. Previous findings suggest that U (which is mainly an 



atomic property) ranges from 1 to 5 eV for TM atoms [39|, |4Cj . For U > 2 eV our results are 
qualitatively similar with small variations in bond lengths and magnetic moments (< 5%). 
Hence, we will discuss data for U = 4 eV in the following, unless specifically mentioned. A 
large U of 8 eV leads to a reduced cohesive energy and to overestimated magnetic moments, 
especially in the case of Mn. All calculations are performed with a plane wave cutoff energy 
of 408 eV. We use a Monkhorst-Pack [4l| 
structure and density of states (DOS) [18|, 



16x16x4 k-mesh for the calculation of the band 

We find that a 3x3 supercell is large enough to avoid drawbacks of the periodic boundary 
conditions. Our supercell contains 36 C atoms and one dopant in all cases, with a unit cell 
length of 7.45 A in the a6-plane. The length of the c-axis is relaxed (~7 A). It was suggested 
that the energy cost associated with the formation of two separated magnetic moments is 
given by the RKKY coupling between two magnetic impurities embedded in a metal, which 
decays slowly with the separation of the impurities in the case of a low-dimensional system. 
Consequently, inclusion of this effect would demand a huge unit cell with several magnetic 
impurities, exceeding reasonable computational resources. We fully relax the positions of all 
the atoms, including the dopants, in order to obtain the minimum energy configuration. We 
continue the optimization until an energy convergence of 10 -7 eV per supercell and a force 
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convergence of 0.04 eV/A is reached. To avoid trapping of TM atoms in local energy minima 
we start have studied different starting configurations where the TM atom is located above 
either the bridge site (center of a C— C bond), or the hollow site (center of a C hexagon) or 
the top site (top of a C atom). The final lowest energy structures are tabulated in Fig. 1. 

III. STRUCTURES AND STABILITY 

Our data suggest that the TM atoms prefer being accommodated in different sites when 
intercalated in aaa and aba hosts. In aaa stacking Fe, Co, Ni, and Cu atoms tend to move 
towards the hollow site, while Mn atoms (smallest atomic number Z) prefer to remain next 
to a C atom and form an almost linear C— Mn— C bond. However, in aba stacking Mn, Fe, 
Co, and Ni atoms favor the hollow site of one C layer, which is the top site of the next C 
layer, see Fig. 1. Cu atoms (highest Z) are slightly off-centered as compared to the other 
atomic species. 

The observed structural variations can be explained by the crystal field splittings which 
apply to the d orbitals of the different TM atoms. If a spatially symmetric negative charge 
is placed around a particular TM atom, the d orbitals remain degenerate. However, being 
itself negatively charged the TM atom experiences repulsion and the orbital energies are 
raised. As in our case the field results from electrons in the C layers on both sides of the TM 
atom, the charge distribution is not spherical and the degeneracy of the d orbitals is lifted. 
The effect is stronger for atoms with partially filled d orbitals (Fe, Co, Ni, and Cu), which 
thus tend to distort away from the trigonal symmetric position and form a tetrahedral-like 
structure with the C atoms. This is clearly visible in the left column of Fig. 1. The average 
C-TM-C bond angles for Mn, Fe, Co, Ni, and Cu turn out to be 180°, 143°, 151°, 144°, 
and 155° in aaa stacking and 143°, 139°, 141°, 141°, and 164° in aba stacking. The Mn atom, 
being exactly half filled, experiences only a weak effect and retains the trigonal symmetry, 
forming a linear C— Mn— C bond. In the aba stacking there is additional negative charge 
due to the C p z orbitals. Here, the shift of neighboring C layers within the a6-plane offers 
the TM atoms a tetrahedral-like geometry. However, the unpaired electron of the Cu atom 
results in some off-centering. 

We address the cohesive energy (E co h) in Table I in order to describe the stability of the 
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different structures. The cohesive energy per atom is defined as 

{E C ell ~ ^ E atom — 

&coh , UJ 

n 

where E ce u is the (spin-polarized) total energy of a cell containing n atoms. E atom and Ed 
are the total energies of isolated C and dopant atoms, respectively. They are calculated by 
placing a single atom in a sufficiently large cubic unit cell of ~11 A lateral length, in order 
to avoid interaction with the periodic image. 

The calculated magnetic moment, smallest distance between the TM atom and its nearest 
neighbor C atom, cohesive energy, change of the in-plane C— C bond lengths, and estimated 
buckling are summarized in Table I. In aaa stacking the TM— C bond length for Mn is 2.00 
A, while it varies for Fe, Co, Ni, and Cu in the ranges 2.03 to 2.26 A, 2.01 to 2.04 A, 1.92 
to 2.18 A, and 1.97 to 1.98 A, respectively. In aba stacking the TM— C bond lengths for 
Mn vary in the range 2.00 to 2.39 A, while for Fe, Co, Ni, and Cu we obtain 2.02 to 2.19 
A, 2.02 to 2.30 A, 1.96 to 2.22 A, and 1.97 to 2.19 A, respectively. Both the bond lengths 
and bonding energies in Table I indicate that the bonding between the TM atoms and the C 
layers is covalent. A significant buckling is observed in the C layers due to the repulsion of 
the TM charge, amounting to 0.185 A, 0.093 A, 0.168 A, 0.107 A, and 0.098 A (aaa stacking) 
and 0.118 A, 0.015 A, 0.079 A, 0.044 A, and 0.098 A (aba stacking) for Mn, Fe, Co, Ni, 
and Cu, respectively. The buckling is calculated as the difference between the z-coordinates 
of the highest and lowest C atoms in a particular layer. The in-plane C— C bond lengths 
range from 1.43 A to 1.45 A in all structures, showing that the sp 2 -hybridized nature of the 
in-plane bonding is retained. 



IV. ELECTRONIC STRUCTURE AND MAGNETISM 

E co h and the bonding strength decrease monotonically with increasing onsite interaction 
U (with a slight exception for Fe), see the top panels of Fig. 2. This has a clear but limited 
effect on the magnetic properties induced by the dopants, as is evident from the bottom 



panels of Fig. 2. We note t 

n 

£/*, consistent with Ref. [18 



iat Mn shows a saturation of the magnetization with increasing 
, whereas for U = 8 eV the magnetization trails off for Fe and 
Co. A growing value of U will affect more prominently the partially filled d systems Fe and 
Co. A choice of U = 4 eV is reasonable for our purpose. 
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In contrast to Mn doped graphene [181 ]. Mn intercalated graphite develops metallic states, 
see Fig. 3. The spin-polarized DOS for [7 = 4 eV shows that the presence of the TM enhances 
the DOS at the Fermi level. Both stacking schemes induce spin-polarization for Mn, Fe, and 
Co intercalation, amounting to 4.10 3.70 and 1.97 \ib (aaa stacking) and 3.80 2.06 
He, and 1.83 \ib (aba stacking). Ni and Cu intercalation does not result in a magnetic state 
for both aaa and aba stacking. For all dopants the induced magnetic moment differs in aaa 
and aba stacking, due to different average C— TM bond lengths. The smaller this length, the 
larger is the magnetization. The spin density induced by Mn, Fe, and Co is addressed in the 
3rd and 5th column of Fig. 1. The maps show localized magnetic moments on the C atoms, 
which reveal characteristic patterns of parallel and antiparallel orientations. When the total 
magnetic moment decreases along the series Mn-Fe-Co, these localized moments reduce in 
magnitude. However, simultaneously a weak spin density develops in the interstitial region 
between the atomic sites which shows a distinct circular wavy shape. This delocalized spin- 
polarization even dominates in the case of Co intercalation. 

A Lowdin charge analysis shows that the Mn atom loses almost 0.5 electrons when it is 
intercalated, with d orbital occupations of ^5.5 electrons in aaa and ^5.6 electrons in aba 
stacking. Due to the trigonal-like crystal field around the Mn atom, the fivefold degeneracy 
of the d orbitals is lifted and the d 3z 2_ r 2 orbital becomes almost fully occupied (in a com- 
paratively weak crystal field). Hence, the 5.5 electrons in aaa stacking yield local magnetic 
moments of 3.5 \±b to 4 /ig, while in aba stacking the 5.6 electrons yield a slightly smaller 
moment. We find that Fe loses almost 0.4 electrons in aaa and 0.3 electrons in aba stacking. 
The occupations of the d bands are 4.8 and 4.7, respectively. A stronger crystal field results 
in partial occupations of the rf 3z 2_ r 2, d XZ) and d yz orbitals and total magnetic moments of 
3.70 \ib and 3.06 \ib for aaa and aba stacking. A very strong crystal field splitting for Cu 
and Ni results in partial filling of all d orbitals and a prominent out-of-plane bonding of the 
d orbitals with the C p z orbitals, resulting in zero magnetic moment. 

Analysis of the DOS reveals that for Mn intercalation and aba stacking the states below 
the Fermi level reveal a strong contribution of the C p orbitals, while the states just above 
the Fermi energy are dominated by the Mn d orbitals. For both aaa and aba stacking, the 
DOS analysis suggests that the Mn atom is in a Mn 2 ~ 6 state, where 5 is a small positive 
number. The partial s, and d occupations reflect the non-ionic character of the bonding. 
In the case of Fe intercalation, for both aaa and aba stacking, the states below the Fermi 
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level are entirely due to the C p orbitals, while the lowest unoccupied states are a mixture 
of C p and Fe d orbitals. In case of Ni and Cu intercalation, we find that the d orbitals do 
not contribute to the electronic states around the Fermi level, which mainly originate from 
C p and TM s orbitals. This explains the non-magnetic behavior of Ni and Cu. In general a 
mixing of C p and TM d states points to a partially covalent character of the TM— C bond. 



V. CONCLUSION 



In view of the possibility to realize a Kondo system by chemisorption of metal atoms on 



graphene 



431 ]. we have addressed the behavior of the TM atoms Mn, Fe, Co, Ni, and Cu 



when they are intercalated in a 3D graphitic network with aaa and aba stacking. While TM 
adatoms on pristine graphene are reported to have bonding energies of some 0.2 eV to 1.5 



eV 



44 



451 ]. we obtain bonding energies between 1.10 eV and 4.23 eV for the systems under 



investigation. The electronic structure of the graphitic systems is modified significantly after 
intercalation due to hybridization of the C p z orbitals with the TM d orbitals. It turns out 
that Mn, Fe, and Co induce spin-polarization in both stacking configurations, whereas Ni and 
Cu result in metallic systems with zero magnetic moment. While graphene is characterized 
by a high mobility of incorporated TM atoms, the strong preferential bonding of intercalated 
TM atoms in our case paves the way to graphitic Kondo systems. 
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FIG. 1: Relaxed structures of the modified graphitic systems with intercalated TM atoms in aaa 
and aba stacking. The induced spin densities are compared to each other. 
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Stacking 


Dopant 


Magnetic moment (/is) 


Itm (A) 


Ecoh (eV) 


E b (eV) 


h (A) 


Ib (A) 




Mn 


4.10 


2.00 


8.93 


4.22 


1.43-1.45 


0.18 




Fe 


3.70 


2.00 


8.90 


3.69 


1.43-1.44 


0.09 


ctctct 


Co 


1.97 


2.02 


8.92 


4.03 


1.43-1.44 


0.10 




Ni 


0.00 


1.99 


8.91 


4.12 


1.43-1.44 


0.17 




Cu 


0.00 


1.98 


8.82 


1.13 


1.43-1.45 


0.01 




Mn 


3.80 


2.19 


8.91 


4.29 


1.43-1.45 


0.12 




Fe 


2.06 


2.17 


8.89 


3.42 


1.43-1.45 


0.01 


aba 


Co 


1.83 


2.26 


8.92 


4.15 


1.43-1.45 


0.08 




Ni 


0.00 


2.15 


8.90 


4.02 


1.43-1.44 


0.04 




Cu 


0.00 


2.04 


8.82 


0.82 


1.43-1.45 


0.09 



TABLE I: Magnetic moment (per unit cell), average distance between dopant and nearest neighbour 
C atoms (Itm), cohesive energy (E co h), in-plane C— C bond length (lc), bonding energy and 
buckling of the graphene layers (Ib)- 
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FIG. 2: Variation of the cohesive energy (E co h) and total magnetic moment (per unit cell) with 
the onsite energy (U) in aaa and aba stacking. 
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E-E F (eV) E-E F (eV) 

FIG. 3: Full lines: Spin majority and minority DOS obtained for the configurations (a) Mn(aaa), (b) 
Mn(aba), (c) Fe(aaa), (d) Fe(aba), (e) Co(aaa), (f) Co(aba), (g) Ni(aaa), (h) Ni(aba), (i) Cu(aaa), 
and (j) Cu(aba), where U = 4 eV. Dashed lines: Corresponding DOS for pristine graphite. 
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